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ABSTRACT 

The  exact  formulas  for  the  bias  and  variance  of  the  product  estimator 
under  simple  random  sampling  are  given .  Its  asymptotic  normality  is 
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rigorously  established  under  weak  and  interpretable  regularity  conditions  on 
the  finite  populations*  No  superpopulation  model  is  assumed.  Hajek's 
projection  method  is  the  key  tool  in  our  proofs . 
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SIGNIFICANCE  AND  EXPLANATION 


Suppose  we  want  to  estimate  the  average  acres  of  corn  per  farm  (denoted 
Y)  in  the  state  of  Wisconsin  based  on  a  simple  random  sample  of  X  values  and 
available  information  on  the  size  of  each  farm  (denoted  X).  In  this 
particular  instance,  X  and  Y  are  positively  correlated.  It  is  common  to 
use  the  ratio  estimator  defined  as  (y-sample  mean)  •  ( x-population  mean)/(x- 
sample  mean)  as  the  estimator  of  y-population  mean.  The  ratio  estimator  is 
very  simple  to  compute  and  is  efficient  for  a  general  class  of  populations. 

In  case  X  and  Y  are  negatively  correlated,  the  product  estimator  defined 
as  (y-sample  mean)  •  (x-sample  mean)/ (x-population  mean)  is  more  efficient  for 
estimating  the  y-population  mean  than  the  ratio  estimator .  Its  computation  is 
as  simple  as  the  ratio  estimator.  Sometimes  we  want  to  find  a  confidence 
interval  for  the  y-population  mean.  Since  no  exact  finite-sample  distribution 
theory  is  available,  we  study  the  asymptotic  normality  of  the  product 
estimator.  Our  regularity  conditions  are  directly  on  the  finite 
populations.  No  artificial  superpopulation  model  is  assumed.  The  conditions 
are  weak  and  interpretable  in  practical  terms  so  that  the  practitioners  may 
find  them  useful. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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THE  ASYMPTOTIC  DISTRIBUTION  OF  THE  PRODUCT  ESTIMATOR 


C.-F.  Wu1'2  and  D.-S.  Chang*'1'3 

1 .  Introduction 

A 

In  sample  surveys  the  ratio  estimator  y  -  yX/x  using  a  suitable 

K 

auxiliary  variate  X  is  one  of  the  most  commonly  adopted  estimates  of  the 
population  mean  of  character  Y  under  study .  In  case  both  variates  Y  and 
x  are  positive  (or  one  positive  and  the  other  negative)  and  the  corretlation 
coefficient  P  between  Y  and  X  has  high  negative  value  (or  high  positive 

A 

value)/  the  product  estimator  Y  *  yx/X  is  preferred.  Throughout  the  paper 

P 

we  assume  the  sample  is  obtained  from  simple  random  sampling  (s.r.s).  It  is 
known  (Murthy/  1964)  that  the  product  estimator  has  smaller  mean  square  error 
than  the  unbiased  estimator  y  ,  if  p  is  smaller  than  minus  one-half  of  the 
ratio  of  the  coefficients  of  variation  of  X  and  Y,  Though  the  asymptotic 
behavior  of  ratio  and  regression  estimators  was  discussed  in  Soott  and  Wu 
(1981),  the  corresponding  results  for  the  product  estimator  were  not 
satisfactorily  given  in  the  literature.  Srivastava's  result  (1966)  is  not 
rigorous,  is  for  infinite  population,  and  regularity  conditions  were  not 
spelled  out.  Had  the  conditions  been  given  there,  they  would  be  unnecessarily 
strong  because  his  approach  is  not  powerful  enough.  The  technique  of  proof 
for  the  ratio  estimator  in  Soott  and  wu  (1981)  does  not  work  for  the  product 
estimator.  The  key  technical  problem  is  that  both  y  and  x  involve  the 
same  random  sample.  As  explained  in  more  detail  in  section  3,  standard 
probability  results  can  not  be  directly  applied  to  this  problem.  The  key 
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technical  tool  employed  in  this  paper  is  Hajek's  projection  method  which  has 
proved  successful  for  many  other  problems  in  statistics .  We  project  the 
product  variable  y  x  onto  a  linear  space  of  random  variables,  thus  obtaining 
a  variable  of  simpler  form.  We  then  prove  the  asymptotic  equivalence  of  the 
two  variables.  The  asymptotic  normality  of  the  latter  variable  was  proved  by 
using  the  Lindeberg-Ha jek  Central  Limit  Theorem.  The  required  regularity 
conditions  are  weak  and  easily  interpretable. 

In  the  next  section  we  derive  exact  formulas  for  the  bias  and  variance  of 
the  product  estimator.  The  asymptotic  normality  of  the  product  estimator  is 
proved  in  section  3,  using  the  projection  method  outlined  in  the  previous 
paragraph.  The  asymptotic  consistency  of  a  variance  estimator  is  established 
in  section  4.  Finally  the  two  results  are  combined  to  give  the  main  Theorem 
3 ,  which  is  followed  by  some  remarks  on  regularity  conditions. 
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2 .  Bias  and  Variance  of  the  Product  Estimator 

Before  discussing  the  asymptotic  results  of  the  product  estimator,  we 
first  derive  the  exact  formulas  of  the  bias  and  variance  of  the  product 
estimator  and  use  the  latter  to  obtain  an  approximate  variance  in  the  large 
sample  case.  The  technique  is  the  same  as  in  Goodman  (1960)  or  Sukhatme  and 
Sukhatme  (1970,  p.  190-192). 

In  population  U  =  {i,2,...,n),  let  f  ■  n/N  be  the  sampling  fraction, 
y  and  x  be  the  sample  means  of  simple  random  sample  of  size  n  of  the 
variates  Y  and  X  respectively,  and  Y  and  X  be  the  corresponding 
population  means .  Denote 

Sc 

and  (1) 

Ca0  =  W*0*0  '  =  °'1'2  * 

Then,  we  have  the  following 
Lemma  1  . 

(i)  E(yx)  =  YX(1  +  —  C,,) 

n  ii 


N 


(N-1)"1  l  (Y  -Y)®(X  -X)B 
i-1 


3O0 


and 

(ii)  Var 


w  {c20%2+2Cii  +  j[2‘i£f>  • 


(n-1 ) (N-1 ) (N-n-1 )  2 

n  (N-2  )  (N-3 )  '  u2t'20^i  1 ' 


(2) 


N-n  2  ~]  1  .  M  +N-6nN+6n  \ 

N  C11  2  (N-2)  (N-3)  22/  * 

—  n  J 


Proof . 


Let 

5x  =  (x-X)/X  ,  Sy  -  (y-Y)/Y  . 
Then,  using  Theorem  2.3  of  Cochran  (1977),  we  have 


(3) 


4 


-3- 


E(yx)  =  YX[1+E(6y6x)J  =  YX ( 1  +  C,4  )  . 

n  11 

Using  results  in  Sukhatme  and  Sukhatme  (1970,  p.  190-192)  and  after  some 
simplication,  we  have 


Var(yx)  =  E(yx)2  -  (E(yx)]2 


(YX) 2E ( (fiy ) 2 (fix) 2+2 (fiy )25x+2fiy (fix) 2+(Sy ) 2+ (fix) 2 


+4fiyfix+2fiy+2fix+1 ]  -  [E(yx)]2 


l™2(-2rl  {c2o%2+Jcii  *1  [2<£ 


-2n ,  , 

T)(C21K12> 


+  (n-lUN-ilUN-n-ll  {c  c  +2c2 

n (N-2 ) (N-3 )  '20  02  1l' 


Nzn  2 

N  11 


2  2 

Jj _  m  N  +N-6nN+6n  1 

+  2  *  (N-2) (N-3)  C22  [* 

n  > 


Applying  the  results  in  Lemma  1,  we  have  the  following  formulas  for  the 

product  estimator  Y  =  yx/X  . 

P 


(i)  E(Y  )  -  Y  =  -  C..Y 

p  nil 


(ill  V4t«p)  -  {C204C02+C,, 


.  1  fNl2n  (n-1  )  (N-1  )  (N-n-  1  ) 

n  jj  N-2  C21^12  n  (N-2 )  (N-3 ) 


,  „  2  ,  N-n  2  “I 

(C20C02+2C11)  -  IT  C11 

2  2 
J_  N  +N-6nN-H6n 

2  (N-2) (N-3)  C22  ’ 
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The  formula  Var(yx)  in  GooOnan  (1960,  (21),  p.  712)  differs  from  oura 

in  that  Goodman' s  treatment  of  product  is  for  sampling  from  infinite 

2 

population.  Me  should  point  out  that  the  coefficient  of  in  his  formula 

(21)  should  be  n“2  instead  of  n"^,  a  minor  error  whidi  disappears  as 


n  ♦  °» 


Following  the  usual  formulation  of  the  central  limit  theorem,  we  embed 


our  finite  population  in  a  sequence  of  populations,  (Uv)  ,  indexed  by  v 

where  ny  and  both  increase  without  bound  as  v  +  *.  And  let 

(sv)  be  a  sequence  of  simple  random  samples  with  sample  size  nv  and  sample 

means  y^,  of  Y  and  X  respectively  in  sv.  Similarly,  let 

-  -  2 

Yy  and  Xy  be  the  population  means  of  Y  and  X,  s  the  population 

variance  of  Y  in  Uv,  and  fy  =*  nv/Nv  the  sampling  fraction.  All  our 

results  are  based  on  the  following  two  standard  theorems  in  finite  populations 

L  P 

(see  Soott  and  wu,  1981).  We  use  - >  and  ■  -  >  to  denote  convergence  in 

distribution  and  in  probability  respectively. 

Theorem  A. 

Suppose  n  +  ®  and  ~  n  +  “  as  v  ♦  "  Then,  under  simple  random 
sampling. 


•'nT  (y  -  Y  ) 
v  v  v  L 


/  i-f  s 

v  vy 


•>  N(0,1)  as  v  +  * 


if  and  only  if  {yvj}v  ^  satisfies  the  Lindeberg-Ha jek  condition 


I 


lim 

T  (6)  (N  -1  )S* 
v  v  vy 


0  for  any  c  >  0  , 


where  Tv(e)  is  the  set  of  units  in  Uv  for  which 


S 

vy 


>  e/JT*  . 

v 


(A) 


'■fW 


Theorem  B. 

Suppose  (y^y)  satisfies  the  condition 

( 1 — f  )S^  /n  +  0  as  v  * 
v  vy  v 

Then,  under  simple  random  sampling, 

(i)  E  |y  -  Y  |  ♦  0 
v  v 


(ii)  y  -  Y 
v  v 


->  0  as  v  + 


The  main  technical  problem  stems  from  the  fact  that  both  y  and  x 

involve  the  same  set  of  indicator  variables  1^,  further  complicated  by  the 

dependency  of  among  i  =  1,2,...,N.  As  far  as  we  know,  no  existing 

results  in  probability  theory  can  be  directly  applied  to  our  problem.  We 

2-  - 

project  the  variable  5^  =  onto  the  space  of  linear  combinations  of 

and  establish  the  asymptotic  equivalence  of  5^  and  its 
projection.  Since  the  projection  variable  is  a  weighted  sum  of  Tif  its 
asymptotic  normality  can  be  established  relatively  easily* 


For  simplicity  we  drop  the  subscript  v  in  all  the  derivations  and 


3 roofs . 


Let  Ii  be  the  indicator  of  unit  i  in  simple  random  sample  s  without 
replacement  of  size  n.  We  can  write  the  centered  version  of  the  product 


2  — 

variable  5  =  n  xy  as 


N  N 


• 5  • •  A  A  xiWi  -  "’iV1  • 


i»i  j-i 


project  h  onto  the  space  Z  of  linear  oombinations  of  {i  }  ^ ,  i.e, 


(z|z  -  £  k  I  ,  k  is  real-valued  constant  for  i  »  1,...,n) 

i-1  1  1  1 


Let  the  projection  variable  be  W.  Then  we  have  (Lehmann,  1975,  p.  362), 


>4 


(6) 


w  -  1  E(n|i  )  . 

1-1 

By  taking  conditional  expectation,  we  have  the  following  properties, 
where  f  -  n/N  and  f1  -  (n-1)/(N-1), 

(i)  Ei^)  -  f  ,  (ii)  "  ffi 

(7) 

(iii)  E(Ij|li)  -  ^-(n-I^,  (iv)  E(IkIj,Ii)  =  f1  N^2(n”2Ii) 

for  i  *  j  *  k  =  1,2,...,N. 

From  (7)  we  have  (for  details,  see  Chang,  1981) 

E"llIl>  ■[,1MS|hAl  +  ttMHMi'Vi  ‘  ■  l  xkV](Ii-f>  ■ 


for  i  -  1,2,...,N  and  from  (6),  we  obtain 


N  N 

W  =  l  a  (I  -f)  -  l  <b  +c  ><I  -f)  (8) 

i-1  11  i=*1 


-8- 


From  k  »  ^  *c,  *  0  and  (9).  we  have 

i*1  i  i 

y  bib3  •  [f-6,is] 

y  ci=i  ■  fow-a)g]i(i,-') 

E 


i > 

i“1 
N  , 

I  < 

i=1 


(N-1)(C20"C02+2C11)  ' 


x  |c  +  (t]-1)(N-2)  c2 
X  |C22  2  11 
N 


*  21=2, «, 2^20^02  * 


no) 


and 


i=*1 


b  c 
1 1 


I  I  ^C.  =  ~2-^-  X2Y2(C^-tC.<-fCOA-K:n^-t-2C,,) 

i*j  1  3  1  (N-2  y 


'12  21  20  02  1 1 


Derivation  of  the  identities  in  (10)  follows  easily  from  definition  (9)  and 
formulas  in  Sukhatme  and  Sukhatme  (1970,  p.  190-192).  Thus,  finally  we  obtain 


Var(W)  -  n3(1-f)X2Y2  ^ 


Lln-UNj_]2 

[n(N-1)(N-2) J  1 


_+2c,,  ) 

20  02  11 


+  2 
n 


(n-1 ) (N-2n)N~ 


2  2  (C2l"C12"C204C02+2C11) 

n(N-1)  (N-2) 


1 

N  (N-2n) 

2 

_  .  (N-1) (N-2) 

2 

n 

(N-1) (N-2) 

C22  2 

N 

11 


+  2^C21+C12+C20-K:02+IFici1+l) 


(11) 


Further,  by  (5)  and  (6),  we  can  decompose 

h  «  W  +  Q 


and  verify  the  identity 


-9- 


Var  (*t)  -  Var(W)  +  Var(Q) 


(12) 


Derivation  of  (12)  is  similar  to  that  in  Lehmann  (1971,  A.  168,  p.  363).  For 
details  see  Chang  (1981,  Lemma  4.2). 

Next,  we  want  to  demonstrate  the  asymptotic  equivalence  between  n  and 


W  as 

follows . 

Lemma 

_2. 

Var (nv)/Var(Wv) 

♦  1  as  v  +  “  provided 

(i) 

n  and  N 

V  V 

both  ♦  *  as  v  ♦  *, 

(ii) 

X  ,  Y  and 

V  V 

C  .  „  a  n  .  are  all  uniformly  bounded 

vap ,  a, p=0 ,1,2 

in  v 

and 

(iii) 

(c  +C  „  +2C  is  bounded  away  from  zero  uniformly 

v20  v02  vll  v 

in  v,  where 

Cvaf5 

is  defined  in 

(1)  for  population  Uy . 

Proof . 

Under  conditions  (i) ,  (ii)  and  (iii),  both  Var(n)  and  Var(W)  in  (2) 

3-2-2 

and  (11)  are  dominated  by  n  Y  X  <c2o4C02+2C1 1 )  as  v  *  “•  This  Proves 

2  . 

Similarly,  under  conditions  (i),  (ii)  and  (iii)  of  Lemma  2,  both 
Var(I^_1  b^d^-f)]  and  Varftf)  are  dominated  by  n3Y2X2(C2Q-*C02+2C11)  as 
v  ♦  °®.  Therefore,  we  have 
Lemma  3 • 

Under  the  same  conditions  as  in  Lemma  2, 

N 

v 

Var(W  )/Var[  £  b  . (I  -f  )]  ♦  1  as  v  ♦  *. 
v  i_1  vi  vi  v 


and 


(13) 


W*  =  t  (Sx  +6y  )(I  -f)  »  j>  (fix  +5y  ) 
i=*1  ill  ies 


f 

X  Y 

* 

T(e)  =  - 

i  : 

—  +  ^  -2 

X  Y 

>  e/var ( W* ) 

• 

Then,  by  applying  the  Lindeberg-Hajek  Theorem  A  to  W*  directly,  we  have 
Lemma  4 . 

Let  and  Tv(£)  be  defined  in  (13)  for  population  Uv.  Suppose  that 

n  ♦  00  and  N-n  ♦  “  as  v  *•  °°.  Then 
v  v  v 


W£//Var(W£)  — — >  N(0,1)  as  v  >  “ 


if  and  only  if  {x  ./X  +  Y  ./Y  }  satisfies  the  Lindeberg-Hajek  condition 

v]  v  V]  v  v,3  3  J 

(A) . 

Combining  Lemmas  2,  3  and  4,  we  establish  the  asymptotic  distribution  of 
the  product  estimator  as  follows 
Theorem  1  . 

Under  s.r.s. 


n  (Y  -Y  ) 
v  vp  p 


->  N  ( 0 , 1 )  as  v  ♦ 


/ 1-f  Y  /C  „.+C  „.+2C 

v  v  v20  v02  vl 1 


provided 


(i) 
(ii) 
( iii) 
and 
(iv) 


n  ♦  •  and  N-n  ♦  "  as  v  +  °°, 
v  v  v 

(x  ./X  +  Y  ,/Y  }  ,  satisfies  the  Lindeberg-Hajek  condition  (A) , 

v j  v  vj'  v  v,j  J 

X  ,  Y  and  C0  on.-  are  all  bounded  uniformly  in  v 
v  v  va$,  a, 0=0,1, 2  1 


(c  „  +2C  is  bounded  away  from  zero  uniformly  in  v. 

v20  v02  vl  1  v 
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Proof 


Decompose 


/n(Y  -Y) 
P 


✓5 


(Y  -BY  ) 
P  P 


'l-f  y/c~*c77+2c7T 

20  02  11 


'1-f  Y  •'C-.+C.  +2C  . 

20  02  11 


n-f  ii 

”  /C,.+C.,+2C,. 

20  02  11 


The  second  term  converges  to  zero  by  condition  (i),  (iii),  and  (iv)  . 
first  term  can  be  rewritten  as 


_ W  +  Q _  =  W+Q _ /var(W) _ 

■^7777 5  /cj««o2t2ci,  ,77v77 s  /s7%75577 


where,  by  conditions  in  (i),  (iii),  (iv).  Lemma  2  and  (12), 

/Var  ( W)  /~3  XY  /C  +C  +2C  ♦  1 

v  '  n  d-f )  20  02  11 

and 

Q//var  (w)  — — — : >  0  as  v  ♦ 

Write  B  =  ^_1  bi(IJ,-f),  C  =  I^1  c^Ij-f)  and  by  (8),  we  have 


W  =  B  #  yvar(B)  +  C 
/Var(W)  /Var(B)  Var(W)  /Var(W) 

Then,  similarly,  by  conditions  (i),  (iii),  (iv).  Lemma  3  and  (10), 

C//Var(W)  — ->  0 

and 


The 


-12- 


/Var(B)//Var(W)  - >  1  as  v  ♦  •  . 

Since,  by  (9)  and  (13), 

B  _  W 
/Var(B)  /Var ( W* ) 

which  converges  in  law  to  11(0,1)  by  Lemma  4.  The  proof  is  completed 
applying  Slutsky's  Theorem. 


-4 

E-  4 


•-I 


i 


t 


( 


D 


4.  Variance  Estimation 


-2 


In  practice,  of  course,  Y  ^C2o+cq2+^C1  1  ^  *8  unknown  an(J  estimated  by 


-2 


its  sample  analogue  y  (c20+cQ2+2c1 1 ) .  Let 


vV  "T1  i'2,c20*°0J*Sot1) 


be  an  approximate  variance  estimate  of  Var(Y  ). 

P 


Then,  we  prove  the 


consistency  of  v(Y^)  as  follows. 


Theorem  2 . 


Under  s.r.s.. 


1-f 


V‘V/('^)’'VC»20,Cv02*2Cv1.>  — >  1  “  ’  * 


vp 


provided 

(i) 

(ii) 


n  +  *  and  N  ♦  ®  as  v 
v  v 


"VfV^Wv.j  “d  tlyv)-Jv)2/Sv201v,j 

condition  (B)  of  Theorem  B. 


(iii) 

and 

(iv) 


X  ,  Y  and  C„a  _  are  all  bounded  uniformly  in  v 

v  v  vaB  ,  a, (5=0, 1,2  1 


2 

{p  }  and  (c  _„+2C  both  are  bounded  away  from  zero 

v  v20  v02  vll  1 


uniformly  in  v,  where  is  the  population  correlation 


coefficient  of  X  and  Y  in  Uv. 


Proof . 


Write 


1 


v(v-p)/(J^)v-2(C20<Co2«c„>-i  - 


'20  02  11 


©'.2^.2 


“  -  2 

~  -  2 

« 

©  c20-C20 

Y 

+  2 

Y 

• 
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I 


•■I 


!# 


where , 


-  2 

-  -  2 

i 

©  C02~C02 

=  C 

02 

Y 

02  x  Y 

®  ^"So  =  C20^S2q 


and 


11  X  Y 


By  condition  (iii),  the  uniform  boundedness  of  Cv20  »nd  Cv02  imply 

that  the  variances  of  x/X  and  y/Y  oonverge  to  0.  Together  with  the 

unbiasedness  of  x  and  y  , 

X  P 


>  t  i  . 

x  y 

Further,  by  condition  (ii)  and  Theorem  B,  we  have 

P  .  1  ^20  _P_ 

>  1  ,  >  1  . 

S02  20 


Finally,  since 


^  -  1  =  ih^=  -  P)  . 


S 


11 


20  02 

2 

which  converges  bo  zero  in  probability  by  the  uniform  boundedness  of  (p  ), 
condition  (ii)  and  the  proof  of  Theorem  4  of  Soott  and  WU  (1981,  p.  101,) 
Hence  we  complete  the  proof . 
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.  k  V '  , 

✓  --•Sax*... 


5 .  The  Main  Result  and  Some  Remarks 


By  combining  theorems  1  and  2,  we  state  the  main  result  of  the  paper. 
Theorem  3. 

Under  s.r.s. 


VrT (Y  -Y  ) 
v  vp  v 


v  /c  +c  +2c 
Jv  v20  v02  vl  1 


— ^>  N<0,1) 


as  v  ♦  *  provided 

(i)  n  ♦  •  and  N  -n  +  *  as  v  ♦ 

v  v  v 

(ii)  {Xv^/Xv  +  Yv^/?v)v  ^  satisfies  the  Lindeberg-Ha jek  condition  (A), 

(iii)  {(Xvj"V2/Sv02}vfj  and  {<Yvj‘?v)2/Sv20}v,j  toth  8ati8fy  condition 
(B)  of  Theorem  B. 

<iv)  Xv,  Yv  and  C^^  Q  1  2  are  a11  tounded  uniformly  in  v 

and 

(v)  {p^}  and  ^cv204Cv02+2Cv11^  toth  are  bounded  awaY  from  zero 

uniformly  in  v. 

Remarks . 

Condition  (ii)  says  that  the  attribution  to  the  total  sum  of  squares  of 
the  (X^/X  +  Y j/Y) 's  from  gross  outliers  about  their  mean  value  should  be 
relatively  small.  Condition  (iii)  on  Y j  and  X^  suggests,  respectively, 
that  the  coefficient  of  variation  of  y  and  x  should  be  reasonably  small. 
Condition  (iv)  is  weaker  than  the  more  typical  assumption  in  the  literature 
(e.g.,  David  and  Sukhatme,  1974;  Krewski,  1978;  Krewski  and  Rao,  1961): 

V  V  cva8  *  X  '  Y  *  CaB  as  v  *  "*  Condition  (v)  is  assumed  to  avoid 
trivial  cases  and  is  satisfied  in  most  practical  situations. 
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